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Abstract 

Inverse scattering transform method of the heat equation is developed for a special 
subclass of potentials nondecaying at space infinity — perturbations of the one-soliton 
potential by means of decaying two-dimensional functions. Extended resolvent, Green's 
functions, and Jost solutions are introduced and their properties are investigated in 
detail. The singularity structure of the spectral data is given and then the Inverse 
problem is formulated in an exact distributional sense. 

1 Introduction 

The equation of the heat conduction, or heat equation for short, 

£$(x)=0, (1.1) 

where the operator 

C{x,d.j:) ^ ~dx2+ dl^~ u{x), X={xi,X2) (1.2) 

for more than 25 years has been known |j to be associated to the Kadomtsev-Petviashvili 
(more precisely, KPII) equation 

{Ut - GuUxi + U^ixixjxi = -3Ux2X2- (1-3) 

The scattering theory for the equation of heat conduction with a real potential u{x) was 
developed in but only the case of potentials rapidly decayi ng a t large distances on the 

x-plane was considered. On the other side, it is well known that ( [l.3[ ) is a (2+l)-dimensional 
generalization of the famous KdV equation: if the function ui(t, xi) obeys KdV then 

u{t, Xi,X2) = Ui{t,Xi + fIX2 - Sfl'^t) (1.4) 



solves (1.5) for an arbitrary constant /i G M. Thus it is natural to consider solutions of (1.3) 



that are not decaying in all directions at space infinity but have 1-dimensional rays with 
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behavior of the type (1.4). The scattering theory for the operator (1.2) with such potentials 
is absent in the literature. Moreover, it is easy to observe that like in the KPI case (see Q) 
the standard integral equation for the Jost solution |^ is meaningless for this situation and 
does not determine the solution itself. In trying to solve this problem for the nonstationary 
Schrodinger operator, associated to the KPI equation, a new general approach to the inverse 
scattering theory was introduced, which was called resolvent approach, see f^-fl^ and 
references therein. In jl^ we developed the scattering theory for the iV-soliton solutions 
given in terms of Backlund transformations of the decaying background potential. These 
results for the simplest case N = 1 are essentially used below. 



Here we apply the resolvent approach to the heat equation (LS) with a potential u{x) 
that is a perturbation of a 1-dimensional potential ui{x) of the kind (L4) by means of a 
potential U2(x) rapidly decaying in all directions 

u{x) — Ui{x) + U2{x). (1.5) 

We introduce and study properties of the resolvent, dressing operators, Jost solutions and 
scattering data and formulate the inverse problem relevant to this case. In fact we consider 



here the simplified version of (L4) in which fi = 0. The generic case is reconstructed 



by means of the Galilean invariance of (1.2). Thus in what follows ui{x) = ui{xi) and, 



moreover, we consider for simplicity the case where ui is the 1-dimensional soliton potential 



(see (3T) below). 

Thus here we apply the Inverse scattering theory to a non-scattering situation since the 
"obstacle" is infinite. Such extension of the Inverse scattering theory results in the new and 
unexpected properties of familiar objects, like the Jost solutions and the spectral data. We 
show that they get specific singularities in the complex domain of the spectral parameter. 
Derivation and description of these singularities are our main results here. The article is 
organized as follows. In Sec. 2 we sketch some general aspects of the resolvent approach that 
are necessary for our construction. In Sec. 3 we present results of embedding the theory of the 
one-dimensional one-soliton potential in two dimensions. Presentation here is based on the 
work |pl|. We describe in detail properties of the extended resolvent and Green's functions of 



operator (|l.2| ) with u{x) — ui{xi). On this basis in S ec. 4 the resolvent of the operator (1.2) 
now with the generic potential u{x) given in (|l.5| ) is introduced and its properties are 
described. The departure from analyticity of the resolvent leads us to definitions of the 
Jost solutions and spectral data and description of their properties (Sec. 5). In this way 
we supply all terms of the Inverse problem with proper meaning in terms of distributions. 
In Conclusion some generalizations and future developments of these results are discussed. 
Main results of this article were announced in our work [|l2|. 



2 Extension of differential operators and resolvent 

In the framework of the resolvent approach we work in the space S' of tempered distributions 
A{x, x' ] q) of the six real variables x = (xi, X2), x' , q e R^. It is convenient to consider q as 
the imaginary part of a two-dimensional complex variable q — Qsjf + iq9 = (^11^2) G 
and to introduce the "shifted" Fourier transform 

A{p; q) = JdxJ dx' e»(P+q»)-'q»-' A(x, x'- q^) (2.1) 

where p £ K^, px = piXi + P2X2 and qs^x = qi^iXi + q23tX2- We consider the distributions 
A{x, x'; q) and A(j)\ q) as kernels in two different representations, the x- and p- representation. 
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respectively, of the operator A{q) {A for short). The compositfon law in the x-representation 
is defined in the standard way, that is 

{AB){x, x']q)^ j dx" A{x, x" \ q) B{x", x' \ q). (2.2) 

Since the kernels are distributions this composition is neither ne cess arily defined for all pairs 
of operators nor associative. In terms of the ji-representation (|]^) this composition law is 
given by a sort of a "shifted" convolution, {AB){p;q) = Jdp'A{p — p';c[ + p')B{p']q). On 
the space of these operators we define the conjugation A*, that in the a;- representations 
reads as 



A*{x,x';q) = A{x,x';q), (2.3) 



where bar denotes complex conjugation, or as A*{p; q) = A{—p; —q) in the p-representation. 
Below we say that the operator A{q) is real if A*{q) = A{q), that in terms of p-representation 
means that A{p; q) = A{—p; — q). The set of differential operators ^{x, dx) — J2 dn{x)d^ is 
embedded in the introduced space of operators by considering the operators D with kernel 
D{x,x') — 'D{x, dx)S{x — x'), where S{x) = d{xi)6{x2) is the two-dimensional (5-function 
and, then, by mapping them in the operators D{q) with kernel 

D{x, x';q) = e^'(^~^')i:)(x, x') = V{x, + q)6{x - x'), (2.4) 

to which we refer as the extended version of the differential operator V. The notion of 
reality for a differential operator D is exactly the condition that its coefficients dn{x) are 
real. 

For the operator ( |l.2| ) the extension L{q) is given by 

L = Lo-U, (2.5) 
where Lq is the extension of C{x, dx) in the case of zero potential, i.e. it has kernels 

Lo{x, x'\q) = [-{dx^ + (72) + {dx, + qi)"^] 5{x - x'), Lq{p\ q) = (iqa - q?)(5(p), (2.6) 
and the multiplication operator U can be called the potential since it has kernel 

U{x,x']q)=u{x)5{x-x'). (2.7) 



Below we always suppose that ii{x) is real, which by (2.3) means that the operator ( [l.2[ ) is 
real also: L* = L. 

The main object of our approach is the extended resolvent M{q) of the operator L{q), 
which is defined as the inverse of the operator L, that is 

LAI = AIL = /, (2.8) 

in the space of operators. Here / is the unity operator, I{x, x'; q) — 5{x — a;'), /(p; q) = 5{p). 
In order to make this inversion uniquely defined we impose the condition that the product 

j dp'M{p ~p';q + s + p')M{p'- q) (2.9) 

exists as distribution in p and q and that it is a continuous function of s in a neighborhood 
of s = when s 7^ 0. 
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Thanks to definitions (2.5), ( |2.6| ), and ( |2.8| ) M is real and in particular, the resolvent Mq 
of the bare operator Lo has in the p-representation kernel Mo(p; q) = 5{p){i(\2 — Qi)^^- As 
function of q it is singular when q — ^(qi), where the special two-component vector 



was introduced. In the x- representation by inverting (2^) we get 



(2.10) 



(2.11) 



For a generic operator A with kernel A{x,x';q) the operation inverse to the extension 
procedure, defined in (2.4) for a differential operator, is given by 



A{x,x';q) 



A{x,x'\q). 



(2.12) 



In contrast with the case of the extended differential operators for which D{x,x';q) = 
D{x, x') = 'D{x,dx)S{x — a;'), in general A{x, x'; q) does depend on q and, moreover, can have 
an exponential growth at space infinity. Therefore A{x, x' \ q) not necessarily belongs to the 
space S' of tempered distributions. The fact that A(x,x'\q) can depend on q will play a 
crucial role in the following. For instance also in the case of the simplest resolvent (2.11) we 
have that the function Mq(x, x' \ q) depends effectively on the variable q and is exponentially 
growing at space infinity. More generally from ( ^.8[ ) we have 

(2.13) 



C{x, dx)M{x, x'\q) = C^{x\ d.^,)M{x, x'; q) = S{x - x'), 



where jC^ is the operator dual to C. The function M(x, x'; q) can be considered a parametric 
(g G K^) family of Green's functions of the operator C. In what follows we use special 



notations for the equalities of the type ( ^.13 ), writing them as 

'tMiq) = M{qfC = /, 



(2.14) 



where C denotes the operator £ applied to the cc- variable of the function M{x, x' \ q) and C, 
denotes the operator dual to £ applied to the x'-variable of the same function. Operation 
( ^.12 ) has no analog in terms of the p- representation. Nevertheless, local properties of the 
kernels in the ^-representation are preserved, and we use the kernels with the hat in what 
follows intensively. 



Thanks to our definitions (2.1) and (2.4) it is easy to see that in terms of the p- 
representation the dependence on the q- variables of the kernels of the extension of a differ- 



ential operator is polynomial (like in the example (2.6)). Correspondingly, the essential role 
in the study of the properties of the resolvent is played by the investigation of its depar- 
ture from analyticity, in particular, by its d-bar derivatives with respect to the q- variables. 
Thus to a generic operator A with kernel A{j)\ q) in the p-representation we associate two 
operators djA with kernels 



(5,^)(p;q) 



g^(p;q) 
9q,- 



J = 1,2, 



(2.15) 



where the derivatives are considered in the sense of distributions. In terms of the objects 
introduced in (2.12) we get by inversion of (|2.lD that 



dA{q) 



2 dq. 



(2.16) 
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Multiplying equalities in (2.8) from the left and right, correspondingly, by Mq we get 
thanks to (2.5) 



M ^Mq + MqUM, 



M = Mr,+ MUMr, 



(2.17) 



Since the resolvent Mq is explicitly given these are integral equations determining the so- 
lution M of (2.8). In the literature (see, say, Q and |^) on the Jost solutions of the heat 
equation some small norm conditions on the potential u are known to guarantee the ex- 
istence of the Jost solutions. So it is natural to assume that under such conditions the 
solution M of the above integral equations exists and is unique (the same for both integral 
equations). In this case the resolvent M can be considered as a small perturbation of th e 
resolvent Mq and this bare resolvent determines the properties of M by means of ( 2.17 ). 
The main problem of construction of the Inverse scattering transform for the operator (1.2) 
is that the potential u{x) in (1.5) does not obey any small norm condition. 

In order to overcome this difficulty we use a so called Inverse scattering transform on 
a non trivial background |l0[ . Let us consider a kind of Hilbert identity, known in the 
standard spectral theory of operators. Precisely, if M [q] is the extended resolvent of the 
operator L{q) with potential u and M'{q) the extended resolvent of the operator L'{q) with 
a different potential u' , then, by ( |2.8| ) we have 

(2.18) 



M' - M = ~M'{L' - L)M. 



Strictly speaking, this follows under the assumption that the product in the r.h.s. is asso- 
ciative. This is a natural assumption since L'{x, x' \ q) — L{x, x'; q) = {u{x) — u{x'))5{x — x') 
and M satisfies condition (2.9). Let now Ci denote the operator ( |1.2| ) in the special case 
where the potential u{x) in (1.^) is purely 1-dimensional, i.e. U2{x) = 0. Let Li denote its 
extension and Mi its resolvent, that is let (cf. (|2.5|)) 



Li — Lq - Ui, 



L^Mi = M^Li = I, 



(2.19) 
(2.20) 



where as in (2/7) Uj{x,x';q) = Uj{x)5{x — x'). Choosing now in ( 2.18 ) L' — Li we get 

M ^ Ml + MiU^M, M = Ml + MU^Mi, (2.21) 



where the second equality is derived in analogy. These equations generalize (2.17) for the 
case where Mi is non trivial and if the resolvent Mi is known they also can be considered 
as defining the resolvent M. If we choose U2 obeying the small norm condition mentioned 
above, we can assume that the solution M of both equations (2.21) exists and is unique. 
Then, thanks to ( |2.20D M obeys (U). Contrary to ( ^.17] ) now M can be considered a 
perturbation of the resolvent Mi. So in the next section we study the properties of the 
resolvent Mi in detail. 



3 One-dimensional potential 

We already mentioned in the Introduction that in this article we deal with the case where 
ui in (^]^) is the 1-dimensional soliton potential 

cosh^ [a {xi — Xq)] ' ^ 
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with a > and xq real constants. In this section we consider the case where the perturbation 
is absent, U2 = 0. We re-formulate in the 2-dimensional space the well known facts about Jost 
solutions of this 1-soliton potential and introduce and study the properties of the resolvent 
and Green's functions in this case. 

The differential equations Ci{x, dx)^i{x, k) — 0, Cf{x,dx)'i'i{x,k) — for the Jost 
solution <I>i(x, fc) and its dual ^'i(x,fc), by using the notation introduced in ( ^.14 ), can be 
re-written as follows 

£'i$i(fc)=0, *i(fc)!Ci=0. (3.2) 

Here and in the following we omit the x-dependence when it is irrelevant. These solutions 
are given explicitly by 

^ k - ^atanh[a{x, - X,)] ^_,,(,),^ 
k — ia 

^ ^fc + zatanhKx,-.o)] ^^^^^ 

k + ia 

where fc G C and the two-component vector i{k) is defined in (^.10|) . They obey the conju- 
gation properties 



= $i(a;,-fc), ^'^(j;, A:) = ^^(a;, -A:) (3.5) 

that are equivalent to the reality condition for the potential ui, and are normalized at 
fc-infinity as follows 

lim e*^('=)^$i(a;,/c) = 1, lim e-'^('=)^*i(a;, fc) = 1. (3.6) 

k — *oo k — ^cjo 

The functions $i(x, /c) and ^i{x,k) are meromorphic in the complex domain of the 
spectral parameter k with poles at k = ia and k = —ia, correspondingly. Thus, these 
functions obey the d-bar equations 

-^^^^tn<P^Jx)5{k-ia), ll^ = in^,^_^{x)d{k + ia), (3.7) 

where we introduced the notations 

*i,aW = res $i(a;,A:), ^^_^{x) ^ ~i res ^{x, k). (3.8) 

k—ia k— — ia 

Explicitly we have 

"'"'"^ ^ cosh [a(x2 — .Tq)] ' ^' ^ cosh[a{xi — Xf^)] ^ ^ 

Let 

c = 2ae^°^° (3.10) 
and ^i^-a{x) and ^'i,a(a^) the values of the Jost solutions in the conjugated points, 

'^i,~a{x)^<^i{x,-ia), ^^^,{x) = ^^{x,ia). (3.11) 
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Then thanks to ( |3.3D , (3.4), and (3.9) the Jost solutions obey in the complex domain of the 
spectral parameter the following scalar products 



j dxi'^^{x,k+p)(^^{x,k) = 27r5(p), 
c [ dx^^^^^^ix)^^ Jx) = 1, 
dxi*i,Jx)$i(a;,/c) = 0, 
dx^ $i__„(x)*i(x,A:) = 0, 
and the completeness relation 

^ / d%$i(a;,fc)«'i(x',/c) + c0(a2_4)$^_Jx)vE'i,Jx') 



pe 



4 < 
4 < a^, 



(3.12) 
(3.13) 
(3.14) 
(3.15) 

^= 5(xi -x'l). (3.16) 



Equations (^j) can be considered as two Inverse problems defining the Jost solution 
and the dual Jost soluti on. The formulation of these problems is closed by giving the 
normalization conditions (|3.6| ) and the following relations: 

(3.17) 



where $i^_a(x) and '^i^aix) are defined in (|3.1l| ). 

These formulae show that the embedding in two dimensions of the Jost solutions of the 1- 
soliton potential is trivial and just mimics the 1-dimensional construction up to the warning 
that, due to their exponential dependence on X2, the functions <i>i_a(x) and '^i^-a{x) are 
not square integrable with respect to the x-variables and, therefore, are not eigenfunctions 
of the operator Ci. 

On the contrary, the 2-dimensional resolvent Mi{q) as well as the Green's function Gi of 
the 2-dimensional operator Ci with the 1-dimensional potential ui are not trivial extensions 
of the corresponding 1-dimensional objects associated to the operator (fc^ -I- — ui(a;i)). 
In terms of the Jost solutions introduced above we can write the kernel of this resolvent 
obtained in ||ll| as 

Mi(x,a;';g) = ^ J dk^[9{ql - - kl^) - 0{x2 - x'^)] ^,{x,k)^,{x' ,k) + 



+ cdia 



ka=qi 
2 „2 



qf) [oia' - q^) ~ e{x^ - 4)] 



(3.18) 



where the hat over the kernel is used in the sense of notation ( 2.12| ). 

Thanks to the equalities ( t3.2| ) and ( 3.16 ) it is easy to check directly that Mi{q) obeys 
the equations 



(3.19) 



that means (cf. (2.14)) that Mi(q) obeys ( ^.20[) and is indeed the inverse of the operator 
Li{q). Moreover, using the explicit formulas ([3.3|), (3.4), and ( |3.9| ) we get that Mi(x, x'\ q) e 
S'{R^), i.e. it belongs to the space of operators under consideration. It can also be proved 
directly that Mi obeys condition (2.9), so it is the extended resolvent according to our defi- 
nition. By means of (|3.5|) we get also that Mi is real, Mj* = Mi according to definition (2.3). 
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We emphasize that in order to prove these results it is not necessary to use the explicit 
formulae for $i and '^i but only their general properties. In fact, if one considers a 1- 
dimcnsional potential ui which, in addition, has a non trivial continuous spectrum one 
gets the same formula for the resolvent Mi. If the discrete (1-dimensional) spectrum of ui 
contains N solitons with parameters aj and Cj (j = 1, 2, . . . , N), then the last term in (3.18) 
must be substituted by the sum of similar terms each corresponding to a value of j. 

Now we describe in detail the properties of Mi{x.x';q). The first term in the r.h.s. of 
) is a continuous function of g = {qi, (72) on the g-plane with discontinuities on the lines 
qi = ±a due to the pole singularities of <l>i(fc) and 5'i(fc). The second term thanks to the 
6 functions has discontinuities on the lines qi = ±a and on the cut q2 — , \qi\ < a. The 
singularities on the lines qi = ±a are exactly compensated among the two terms. Thus 
the kernel Mi(x, x'; q) is a continuous function of q with a discontinuity on the cut 52 — a'^, 
\qi I < a. This discontinuity is specific of the potential ui{x), or more generally of a potential 
with discrete s pectr um and it gives the essential difference of Mi with respect to the bare 
resolvent Mq ( ^.11 ). Let us underline that, in spite of the fact that £1 applied to the 



term with 9{a — 52) that causes this discontinuity gives zero, this term cannot be omitted 



in ( ^.18 ), since only thanks to the fact that 0{x2 — x'2) and 0(a^ — 92) have opposite signs 
the kernel Mi{x, x'; q) = e~i^^~'^ ^Mi{x, x'; q) is a tempered distribution with respect to the 
a;-variables. 



The kernel Mi(p;q) in the p-representation is not an analytic function of q. By ( 2.16 ) 
the d-bar derivatives of Mi with respect to are proportional to dMi/dqj and for the 
latter we get from (3.1g) equalities 



dMi{q) _i 



dqi 

dMijq) 
dq2 



d% M(4cj(fc) - 92) $l(fc) ® ^'l(fc), 



-1 
2^ 



dk^S{£2Q{k) - q2) ^i{k) ®^fi{k), 



(3.20) 
(3.21) 



where ($i(fc) ^ ^'i(fc))(a;, a;') = $i(x, A:)^'i(a;', fc) is the standard direct product and where 
by the above discussion we consider q2 ^ a? ■ For the discontinuity along this line we get 



MM 



-MM 



q2=a^-0 



= -ce{a^ -ql)<^>l,^a®'i'l,a■ 



{3.22) 



We have to study now the behavior of Mi{q) at the end points of the cut, i.e. when 
q ^ (±a, a^). Firstly, it is convenient to subtract to Mi{q) its value, say, on the upper or 
lower edges of the cut: 



9i = lim MM 

ij2=a ±0 



\qi\<a 



(3.23) 



Since ^i{k) and 5'i(fc) are analytic for < a, we deduce from (3.18) that gf are inde- 
pendent also of qi and their kernels equal 



gt{x,x') 



(^2 ^2) 
2tt 



da<i>i{x,a)'i'i{x' ,a) T c9{±{x2 - x'2))^iJx)'i'iJx'), 



(3.24) 



where / da denotes integration along the whole real axis. Now extracting explicitly from 
the first term in the r.h.s. of (3.18|) the contribution coming from the poles of $i(fc) and 
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we get that, say, difference Mi{q) — behaves at points q — (±a, a^) as 



MAq)-gi 



arccot 



+ o(l), g-(±a,a2). (3.25) 



Thus Mi{q) is bounded but discontinuous at q = a^), while its regular part, , is the 
same for both these points. 

Now it is easy to see that the discontinuity of the resolvent along the cut q2 = a^, 
\qi\ < a and the ill definiteness at the points q ~ (±a,a^) are the result of embedding the 
1-dimensional potential in the two dimensional space. Indeed, the resolvent of the Sturm- 
Liouville operator — ui(a;i) — q2 is obtained from Mi{q) by means of the operation 
/ dx2 e-«^(^^-^2)Mi(a;, x'; q). By jsl^ ) and Q, 0, Q we get the standard expression 
for the 1-dimensional Green's function with a pole at q2 = a^. 

We already noted that M^{x,x';q) defines a family of Green's functions. Among them 
we expect should play a special role those obtained considering the values of q belonging to 
the support of the defects of analyticity given in (H^), ( |l2l| ) and in ( |3^ ). We consider, 
therefore, the Green's functions 



Gi{x,x',k) = M^{x,x';q) 
Gt{x,x';k) = M^{x,x';q) 



(3.26) 
(3.27) 



where A; S C is the spectral parameter and we denote qi = /cq, (see (2.10)) in order to meet 
the standard notation. From these definitions it follows directly that 



C iGi(fc) = Gi(fc) £ 1 = /, £ iGf (fc) = (fc) £ 1 = /, 



G,ik) 



fc»=0 



Gi(0), Gf(fc) = G±(zfc3), 



(3.28) 
(3.29) 
(3.30) 



i.e. Gf{ k) are independent on fcjf and then inside the strip they coincide with gf introduced 
in (]3!23|), 



G^" {x^ X 5 Aj) 



: gf{x,x'). 



(3.31) 



As well from (3.18) we get the representations 



Gi{x, x' ,k) = 



= — / dk' 

k' -k 



K\) - eix2 ~ x'2) <fi(x,fc')*i(a:',fc') + 



+ ce{a- \kQ\)9{x'2 ~ X2)^,,.a{^)^U^'), 



(3.32) 



Gt{x,x',k) = 



= - I dk'. 



(fc^)2) - e{x2 - x'2) <i>i(x, k')^,{x', k') T 



T Ceia" - kl)e{±{x2 - 4))$l,-aW«'l,a(^') 



(3.33) 
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The first of these equahties shows that the cut of the resolvent at (72 = a^, |9i| < a is not 
inherited by Gi{k) (in contrast to the case of the nonstationary Schrodinger equation, as 
mentioned in the Introduction) and that Gi(fc) is discontinuous only at the points k — ±ia. 
Its behavior in the neighborhoods of these points follows from ( 3.25| ) and reads as 



G-^{k) — gi < arccot 



*l,-a®*l,a + 0(l), kr.±ia. 



(3.34) 



Also the Green's functions Gf{k) are discontinuous only at feg. = ±a and one gets thanks 
to ( ^3.25 ) that for k ^ ia or k ^ —ia 



Gf(fc) = .gr-c 



1 ±1 



*l,a + 0(l)- 



(3.35) 



Notice that these functions Gf{k) coincide when |/cq| > a and are independent of fccj (and 
then of k by ( |3.30 )) when {k^l < a. On the borders of these strips they have the discontinuity 



Gf(fc) 



outside 



Gf(fc) 



inside 



C 

±-$1 



l,o- 



(3.36) 



Taking into account the discontinuous behavior of the Green's functions we see that equalities 
of the type Gf{ia) = Gi{ia) and Gf{—ia) = Gi{—ia) have no meaning in our case. Thanks 
to ( [3.32 ) and ( |3.35[ ) we have only that 



lim lim G, (k) — g, 

IfcaHa-Ofcjj^O 



(3.37) 



where the limiting procedure must be performed in such way that |/csR|/(a — \k<^\) — > +0. 

In order to complete the study of the Green's functions we mention that Gi{k) obeys 
the standard equalities 



lim (-2ifc)-^ e'^('=)(^-^')Gi(x,a;', k) 



dG^{x,x\k) sgnfcsj 



dk 



2tt 



S{x — x' 
fc)«'i(a;',-fc). 



(3.38) 
(3.39) 



The first of them follows either from the differential equations ( 3.28| ), or from the integral 
representation ( 3.32| ) and properties (3.6). The second one also follows from ( 3.28| ), or it 
can be derived from (3.26) by means of (3.20) and (3.21). This equality must be understood 
in the sen se o f distr ibutions and we see that the discontinuity of Gi(fc) at points k — ±ia 
leads (by ( t3.3[) , (3^)) to the pole singularities of the r.h.s. at these points. In view of ( 3.39 ) 
in what follows we refer to Gi(fc) as the Green's function of the Jost solutions. 



4 Resolvent of the perturbed L-operator 



Now we consider the general case of the operator ( [l.2[ ) with potential given in (|T5|) , where 
1*2 (x) is a real function of two space variables, smooth and rapidly decaying at space infinity. 
The extended resolvent M{q) is determined by (one of) Eqs. ( ^.21 ) and we need to study its 
analyticity properties first. The increment M{p; q + s) — M{p; q) of M can be obtained from 
the Hilbert identity ( 2.18| ) where prime means the increment s of q. We have M' — M = 
-M'{L[ - L^)M and, then, using (|2.20D 



M' - Af = M'L[{M[ - Ml 



(4.1) 
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Thus for the d-bar derivatives with respect to we get 

d^M = {ML^){d^M^){L^M), J = 1, 2, 



(4.2) 



in the region where A/i is continuous, i.e. for q2 ^ c? . In terms of the objects introduced 
in (2.12), we obtain 



d^M = M C ,{d^Mi) C lAf, J = 1, 2, 



(4.3) 



where we used that Li{x, x' \ q) = Cid{x — x') and took into account that when kernels with 
hat are considered the multipUcation by Ci is no more associative and it is necessary to use 
the arrows to indicate the correct order of operations (cf. ( p. 14 )). Now, thanks to (2.16) 
and using (3.20) and (3.21) we get for q2 ^ 



dM{q) _ i 
da, TT 



dM{q) 
dq2 



ka=q 
-1 
2^ 



dk^ M(^23(fc)-<72)(f>(fc)®Vl/(fc), 



fca=?i 

where $(fc) and '^{k) are defined by 

$(fc) = G'(fc)'£ i$i(fc), *(fc) ^ *i(fc)^iG(fc), 

with 

G(x,x\k) ~ M{x,x' ]q) 

9=^3 (fc) 



(4.4) 
(4.5) 

(4.6) 
(4.7) 



More expUcitly, say, the first of equations (4.6) stands for <J>(x, fc) = J dx' (^Cf{x' ,dx')y. 
Gi{x, x' , k))^i{x' , k). The function G{k) with kernel G{x, x' , k) defined in ( [4. 61 ) satisfies the 
differential equations 



CG{k) ^G{k)C =1, 



(4.8) 



which can be obtained as a direct reduction of ( 2.14 ). Therefore, G{k) is a Green's function. 
Since the reduction is the same used in (3.26) for getting Gi{k) from Mi we derive from (2.21) 
that this Green's function obeys the integral equations 



G(fc) = Gi(fc) + Gi(fc)[/2G(fc), G(fc) = Gi(fc) + G(fc)L/2Gi(fc). 



(4.9) 



Again, as in Sec. 3, thanks to (4.7) and (4.4), ( |4.5| ) we get the d-bar derivative of the Green's 
function in the form 



dG{k) _ sgnfcjj 



dk 



27r 



$(-fc) (g) ^'(-A;), 



(4.10) 



where $(fc) and \E'(fc) are defined in (4.6). These objects due to their definition and (4.9) 
obey the integral equations 



$(/c) = $i(/c) + Gi(fc)t/2$(fc), *(fc) = *i(fc) + ^(k) U^G^ik), 



(4.11) 
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where again the first equation more exphcitly reads as k) = $i (.t, k)+J dx'Gi [x, x' , k) x 
U2{x')^{x' , k). It is clear that the differential equations 



C^k) = 0, *(fc) £ = 0, 



(4.12) 



hold and, therefore, we can consider $(a;,fc) and ^'(a;,fc) as the generalization of the Jost 
solutions to the case where the perturbation U2{x) is different from zero. Let us mention 
that thanks to these definitions we succeeded to avoid the indeterminacy in the definition 
of the Jost solutions discussed in the Introduction. Below we study the properties of the 
Green's function and t he Jo st solutions in more details and discuss the singular structure of 
the terms involved in ( 4.10 ). Now let us mention the following standard properties 



2i lim 

k 



im kd^^ fe*^('^)(^~^')G(a;,a;',/c)) =5{x-x'), 

— >oo \ / 



G(fc) = G(-fc) = G(fc), 
^{x, k) = $(a;, -/c). 



^{x, k) = ^{x, -fc), 



(4.13) 

(4.14) 
(4.15) 



that can be obtained by means of the integral equations ( |4.9| ) and properties ( |3. 
and ( |3.38 ) for the Green's function Gi(fc). 

Till now we studied the departure from analyticity of the resolvent in the case 52 ^ a^- 
Since the resolvent Mi{q) is discontinuous along the line 52 = 0? (see (3.22)), the integral 
equations (2.21) suggest that also M{q) has a discontinuity. Let us denote the limiting 
values on the two edges of the line by 



Then from the Hilbert identity ( |4.l| ) we derive that 

M+{q) - M-{q) = M^{q)LM{M+{q) - iq))L,iq)M^ (q), q^ = 



(4.16) 



(4.17) 



where the l.h.s. is independent of the choice of the sign in the r.h.s. In analogy with (3.27) 
we introduce the two Green's functions 



G^{x,x';k) = M{x,x']q) 



(4.18) 



and rewrite (|T7|) in these terms as G+{k) - (k) = {G^{k) C i){GJ (k) - G^{k)) x 
(^iG=F(A:)). Then by ( ^^ and (|^ we get 

G+(fc) - G-(k) = -ce{a^ - kl)<i>^{k) ® ^"^{k), (4.19) 

where the new solutions (cf. ( [4.6D )were introduced: 

<I>±(fc) = G±(A:)tCi<i>i,_„ ^^ik)^^,^,t,G^ik). (4.20) 

Following properties of G^{k) it is easy to show that these Green's functions obey the 
following differential and integral equations and reality condition 



'^G^{k) = G^{k)'c = I, 

G^{k) = Gf (fc) + G^{k)U^G^{k), 



G±ik)^G^{k). 



(4.21) 

G±(fc) = Gfik) + G^{k)U^G^{k), (4.22) 

(4.23) 
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By definition they are independent of fcsjf and by the corresponding properties of G^{k) we 
have that G~^{k) = G~{k) when \kc^\ > a and they are independent of kcj when |fccj| < a. 
By (4.20) and (4.21) we get that ^"^{x, k) and 5'^(fc) are solutions of the heat equation with 
potential (f^), 



/:$±(fc) = 0, ■^^{k)C =0. 



(4.24) 



Integr al eq uations for these solutions follow by applying operations ( 4.20 ) to the equa- 
tions (^4.22|) 



(4.25) 



Let us also mention that thanks to ( 4.23 ) these solutions are real and are in depe ndent of k 
inside the strip < a, due to the corresponding property of G'^(fc) and ( 4.2C ). Since in 
the following we use intensively the Green's functions and these solutions inside the strip it 
is convenient to introduce the following specific notations 



g^{x,x') ^ G^ix,x',k) 



\ks\<a 



and also 



(/)=^(a;) = ^^{x,k) 



I I <^ 



|fcQ|<a 



(4.26) 



(4.27) 



Equality ( 4.19 ) enables us to find relations between solutions ( |4.20|) . Let |fccj| < a, then 
applying, say, C _^ to this equality from the right and using ( |4.20| ) we derive that 



where 



(1 + A)0+ =r, 



(4.28) 



(4.29) 



Explicitly \ = cjdxj dx'^ ^^^{x){C^{x, d^)Cf{x' , d^')g- {x, x')^^^^^{x'). By (|^ this 

constant is real and thanks to ( 4.20| ) it is also equal to A = c(4'i^a ^ 10 ) — c{%Ij^ C i$i^_a)- 
Inserting here £i = £ + C/2 we get by (|4.21|) and ([4.24|) that 



(4.30) 



or A = c(^'i^at/2(/> ) = c(V' C/2$i,-a), whcrc we also used C i$i,_a = and ^i^a >C 1 = 
tha t follo ws from ( |3.2| ) and (3_^). Next, applying to ( 4.28 ) ^'^ ^ £ from the left and again 
by ( |4.20| ) we get 



(l + A)[l-c(vI/i^,£ig+£i$i^_J]^l, 



(4.31) 



where a new constant {^i^a C ig+ Ci<^>i^-a) = («'i,aC/2$i,-a) + (*i,aC/2ff+C/2$i,_a) (cf. 
( 4.3C| )) appeared. Since we chose U2 to be rapidly decaying at infinity all terms must be 
finite. Then 1 + A and, more precisely, taking into account that for U2 ^ also A — » 
we have that 



1 + A > 0. 



(4.32) 
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Summarizing, we get the following relations: 



1 + A' 



+ _ 



l + A' ^ 1 + A' 

ce{a-\h.\) 



G+{k) = G-{k) 



1 + A 



(j) €5 "0 



(4.33) 
(4.34) 
(4.35) 



Here ( [4.33 ) is just (4.31), the first equality in (4.34) is ( 4.29| ) and the second equality is 
deriv ed by analogy, and (4.35) follows from ( 4.19] ) thanks to (4.34). In their turn (4.33) 
and ([4.34I) follow from ([4.351) thanks to (|4.2(]| ) and (|!2| 



We have shown in (3.36 ) that the Green's functions Gf (fc) are discontinuous at fcs 



and — —a. By ( 4.22 ) we deduce that G^(fc) have t he s ame behavior. In order to study 
this discontinuity we use, as above, the Hilbert identity (4.1) where M — AI{q), M' — AI{q'), 



etc. We choose q2 — q'2 — a'^ i: 0, qi — a — e, q[ = a + e, and in the limit e +0 we use the 



hat notation ( 2.12 ) and definitions ( 3.27 ) and ( 4.18 ) of the Green's functions. Then we get 
G±(i(a + 0))-g± = G^{i{a + 0))'Ci{Gt{i{a + 0))-gt)~^ig^, where again (|3^ ) and (|^ 



where used. Now by ( |3.36 ) for the discontinuity of the unperturbed Green's functions we 
obtain 



G=^(i(a + 0)) -g^ = ±-$±(i(a + 0)) ® V"* 



(4.36) 



where notations ( 4.20 ), ( 4.2(: ) and ( 4.27 ) were used. Applying £ from the right 

and ^'i^a C 1 from the left in analogy with the derivation of (4.35) we get by ( 4.2^ ) that 



c(«'i,a£iG=^(i(a + 0))/:i$i _q) = 2A(2 + A)"\ that is finite due to ( [4.32D . Then omitting 
details we derive the equalities G^(i(a + 0)) = G^(— i(a + 0)) = 
that, say, for the bottom sign can also be rewritten in the form 



c(2 + A)-V" 



G-{k) 



c9(\k^^ 



2 + A 



) Tp + 0(1), k ~ ±ia, 



(4.37) 



where we took (4.26) into account. 



5 Properties of the Jost solutions and Inverse problem 



In this section we complete the investigation of the properties of the Jost solutions by 
describing their behavior at the points k = ±za. Formulae ( ^^ ) suggest to study first the 
behavior of the Green's function G(fc). We expect that it is ill defined at these points, so in 
order to describe this behavior we compare G{k) with some well defined Green's function, say, 
g~ . For this aim, as we have already shown, relations of the type (4.19) can be very useful. 
In order to derive them we start again from the Hilbert identity (4.1) where M' = M{q') 
and M = M{q) and we choose q' = £a(fc), qi = k^ 52 = -0 (sec ( |2.10D , (O) a nd (|4.16| )). 
Then, passing to the objects with hats by ( 2.12| ), recalling definitions (4.7), (4.18) and 
keeping only the leading term in the neighborhood of fc ^ ±m we get 



G(fc) - G-{k) = G{k) C i(Gi(fc) - G^{k)) C ^G' (k) + o(l), k - ±ia. 



14 



Inserting the explicit singular behaviors of G'i(fc), (k) and G (fc) at k ~ ±ia given 
in (IsJI), (Osl) and (pTh, wc have 



G(fc) - g- 



2 + A 



(8-0 



- c arccot , ' + ^' - 



G{k)C,^,,^,(E)^-{k)+o{l), 



where in the last multiplier the definition of ^' (fc) in ( 4.20 ) was used. Again by ( 4.20 ) 
and (p7|) VI/- (fc) - — 2+X 
by (|4.27|). Then 



2+A6i(a Ifcal) ^ ^ where as always ip denotes 5* (fc) for |fcQ.| < a 



G{k) - .g- 



2 + A ^ 
2 + A6l(a- Ifccjl) 



c I arccot ,, ' + ^ 

TT 2 



2 + A 



G(fc)/:i$^^_^ ®^-+o(l). 



(5.1) 



Thus in order to get the behavior of G{k) we need to find that of G{k) C i$i^_a, which 
follo ws by applying to ( ^ ) operation £ i$i,-a from the right and using again ( 4.20 ), ( 4.27 ) 
and ([4.29D. Then 



G(fc)/:i$i,_, 



A{k) 



o(l), k ±ia, 



(5.2) 



where we denoted for shortness 



A{k) = TT + A arccot ■ 



(5.3) 



This function is real, positive thanks to (4.32) and discontinuous at A; = ±ia. Now insert- 
ing (5.2) in (5.1) we derive finally that 



G{k)^g- 



arccot 



a - |fccj| 



Aik) 



® ip + o(l), k ^ ±ja. 



(5.4) 



Applying to (5.4) from the left the operation 'i!^ ^ C ^ and recalling the definitions (4.20) 
and (4.27) we derive 



*l,a^lG(fc) 



Trip 

m 



+ o(l), k ±ia. 



and by (4.29) also 



c^.,C,G{k)C,^,^_, 



ttA 

m 



o(l), k ^ ±ia. 



(5.5) 



(5.6) 



Correspondingly, we get for the behavior of the Jost solutions in the neighborhood of k = 
±ia, thanks to (3.3) and (5.2) 



$(fc) 



inca 



A{k){k-ia) 



+ 0(1), kr^ia, 



(5.7) 
(5.8) 



15 



and analogous relations for ^'(fc). 

Now we are ready to consider the d-bar derivative in the sense of distributions of the 



Jost solution, say, $(fc). Let first k ^ ±ia. Then we use ([3. 3D, (|4.6|), and (4.10) to derive 

- <f>(-fc)r(fc), k^±ta, (5.9) 
where the Spectral data are defined as follows 

r(k) = [^,i~k)tMk)^i^iik)). (5.10) 



Thanks to (3.3), (3.4) and ( |5.6D we get the singular behavior of these spectral data in the 
form 

<k) = ^, '4?^^^^^Wt7TT + o{1), k^±ia, (5.11) 

^ ^ 2(fcjp + i|fccj| -ia)A{k) ^ ^' ^ ^ 

i.e. in both points it has a pole singularity multiplied by the discontinuous function A{k). 
Taking into account that the singular behavior of $(— fc) is given by the denominator 
[k — ia)A{k) at point k — ia and by A{k) at point k = —ia we see that the r.h.s. in (5.9) 
is integrable at the latter point but it has a singularity sgnfcjflA: — ia\~'^ A{k)~'^ at point 
k = ia, that is not integrable. On the other side, $(fc) is locally integrable for any k so 
$(a;, fc)e*^'-'^^^ is a Schwartz distribution with respect to k. Thus its d-bar derivative in the 
sense of distributions exists and can be defined in the standard way. Let f{k) be a test 
function that properly decays at infinity (we are not interested in the exponential growth 
due to the multiplier e~*^(*'')^ now). Then the d-bar derivative of $(fc) is defined as 

\k±ia\>e 

where in the last equality we again used the property of local integrability of <i>(fc). Inte- 
grating by parts for e > we can use ( ^.9D and we have 

df{k) 



lim / (Pk $(fc)- - 
e^O J dk 

\kizia\>€ 



lim (h dfc$(fc)-Mim / (fk^{-k)r{k)f{k), 



2i e^O / £^0 

\k—ia\—e \k—ia\>e 

where we omitted the term since thanks to ( ^3.^ ) it gives zero in the limit e ^ 0. 

Thanks to (5J) and ( ^.^ ) both limits in the r.h.s. exist. To be more precise let us introduce 
the distribution 



, 2^^gn%/(fc^^ r ^ sgnfc^/(fc) 

P-"^- ' |fc - za|M(fc)2 - £^™o J ^ > - zaPA(fe)2 " ^^'^^^ 



\k—ia\ >£ 

Notice the presence in the numerator of sgnk^ that guaranties existence of the limit. We 
used the principal value (p. v.) notation in analogy with the one-dimensional case. It can be 
checked directly that 
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where S is some real positive parameter and the second term in ( 5.13| ) is independent on 
the choice of S. In the case where a distribution has singularities of this form at some finite 
number of points oi, 02, etc., we use the same notation for the integral assuming that either 
the cutoff procedure in (5.12) or the subtraction procedure in ( 5.13| ) is performed at each 
point. Of course, the parameters sj and Sj must be chosen in such way that corresponding 
discs do not overlap. 
Let us denote 

\k—ia\—s 

SO that can be considered as an extension of the definition of residuum to the case 
in which th e p ole singularity is multiplied by a function discontinuous at the same point. 
Thanks to (5/7) we get that this limit also exists and equals 

log(l + A) ^_ 



= c- 



A 



Thus, summarizing all above definitions we get that 



dk 



= <^{-k)r{k) + i7r$^(5(fc - ia), 



(5.15) 



(5.16) 



where ^(x, ~k)r{k) is now a distribution in k defined by the p. v. prescription given above. 



By (5.7), ( |5.8| ), and (5.11) it is integrable aX k = —ia, but it behaves as \k — ia\~'^ (k) in 
the neighborhood of the point k = ia. 

Equation (5.16) supplies us with the first equation of the Inverse problem. In order to 
close it we need the analog of the first relation in (3.17), where it is stated that the residuum 
of the function is proportional to its value in the conjugated point. But in our case <I>(fc) is 
discontinuous at point k = —ia, so again some modification of the notion of "value" at this 



point must be given. Following the procedure used in (5.14) we can define it as 

1 . f dk 



lim 

27ri e^o 



■$(fc). 



\k-\~ia\—6 



Thanks to ( |5.q ) this hmit also exists and equals 

log(l + A) 



so that by ( 5.15| ) we have 



A 



= c<I> 



(5.17) 



(5.18) 



(5.19) 



that shows that the parameter c is not modified by the perturbation. This equality closes 
the formulation of the Inverse problem (5.16). Finally, taking into account the asymptotic 
behavior of ^{x,k) and (5.17) and ( 5.19| ) we can formulate the Inverse problem as the 
following system of integral equations: 

$(a;,fc) = e~'^('=)^ + 



I k — k 



My „t(l(ia)-l{k))x 



k 



<i>,(x), (5.20) 



c 

1 



TT J k + ia 2a ° 



p.v. / d^k ■ 
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The integrands in the r.h.s. of these two equations are not locally integrable, respectively, the 
first at k ^ ia and the second at fc = ±ia. Correspondingly, their integrals are regularized 
by means of the principal value prescription, as in ( ^.12 ) or (5.13), at k = ia and at fc = ±ia. 
The potential is reconstructed by means of 



TT J OXl \ 



'k)r{k) 



OXl \ 



(5.22) 



6 Conclusion 



In this article on the basis of the resolvent approach we gave a detailed presentation of 
an extension of the inverse scattering theory for the heat operator to the case where the 



potential (1.5) is a perturbation of the 1-dimensional one soliton potential ui{xi) (B.l) by 



means of a smooth, decaying at infinity function U2{x) of two space variables. To our 
knowledge this is the first time that inverse scattering theory is applied to a non-scattering 
situation, i.e. a situation with an infinite obstacle. As a result of our investigation we 
prove d t hat under such a perturbation the Jost solutions get specific singularities (5.7) 
and (5^) on th e com p lex p lane of the spectral parameter fc. We demonstrated that the 
d-bar problem ( 5.16| ), ( ^.19| ), while looking familiar for a potential whose spectrum has a 
discrete and continuous part, needs a substantially modified approach due to the singularity 
structure of the spectral data given in ( 5.11 ). It was necessary to establish the meaning 
in the sense of distributions of all terms involved in this problem, in order to be able to 



formulate the Inverse problem as the system of integral equations (5.20), (5.21). It is easy 
to check that the singular behavior of the spectral data and Jost solution as given in (5.11) 
and (5.7), (5.8) is compatible with this Inverse problem. On the other side, it is necessary 
to prove that the potential u{x) reconstructed by means of (|^ is of the type (jTI). We 



plan to address this problem in a forthcoming work. 

Another open problem is the application of these results to the KPII equation (1.3) 
itself. In particular, investigation of the time asymptotics of solutions with initial data of 
the type (1.5) must be performed. Let us mention only that the singular behavior ( ^.11 ) of 
the spectral data is preserved under evolution ( |1.3D . Indeed |Q, the time dependence of the 
spectral data is given as 



r{k,t) 



Thus we get that 



a — const, 



A — const 



(6.1) 



(6.2) 



also with respect to time. 

In Sec. 3 we mentioned that the above construction can be easily generalized to the case 
where the potential ui{xi) is a pure iV-soliton 1-dimcnsional potential. At the same time 
our approach also admits straightforward generalization to the case where ui{x) is not a 
function of one space variable but the result of application of the Backlund transformation 
to a generic background 2-dimensional poten tial uo {x) d ecaying on the x-plane. Then the 
inverse problem is again given by Eqs. ( 5.2C ) and ( ^.2l| ), where the spectral data r(fc) are 
replaced with 



r{k) 



(fc + ia){k + ia) 
(fc — ia)(k — ia) 



roik), 



(6.3) 
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where r{k) is of the type ( ^.11 ) and ro(fc) are the spectral data of the potential uo{x) 
(see 0). 

The theory of the heat equation with respect to the nonstationary Schrodinger equation 
is in some respects simpler and in some other respects unexpectedly more difficult. As we 
have shown, under perturbation the Jost solution get singularities more complicated than 
poles, but this solution has no additional cut in the complex domain, in contrast with the 
nonstationary Schrodinger case as discovered in jl^. On the other side the generalization 
of this scheme to the case of multi-ray structure of the potential u{x) meets with essential 
problems, first of all due to the fact that the resolvent (or Green's function) of the heat 
equation even of a 2-soliton (generic) potential is unknown in the literature. This problem 
also needs future development. 
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